Edge critical currents of dense Josephson vortex lattice in layered superconductors 
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We calculate the field dependence of the critical current for the dense Josephson vortex lattice 
which is created by large magnetic field B applied along the layers of atomically layered supercon- 
ductors. In clean samples a finite critical current appears due to the interaction of the lattice with 
the boundaries. The boundary induces an alternating deformation of the lattice decaying inside the 
sample at the typical length, which is larger than the Josephson length and increases proportional 
to the magnetic field. The exact shape of this deformation and the total current flowing along the 
surface are uniquely determined by the position of the lattice in the bulk. The total maximum 
Josephson current has overall 1/B dependence with strong oscillations. In contrast to the well- 
known Fraunhofer dependence, the period of oscillations corresponds to adding one flux quantum 
per two junctions. Due to interaction with the boundaries, the flux-flow voltage for slow lattice 
motion also oscillates with field, in agreement with recent experiments. 

PACS numbers: 74.72.Hs,74.60.Ge 



I. INTRODUCTION 

Atomically layered superconductors, such as 
Bi 2 Sr 2 CaCu20 K (BSCCO), behave as stacks of Joseph- 
son junctions. This so-called intrinsic Josephson effect 
has been a subject of intense research in the past decade 
(see reviewsa'EI). In spite of atomic-size separation 
between the neighboring junctions, in zero magnetic 
field the system approximately behaves like an array 
of independent junctions. A weak dynamic interaction 
between the neighboring junctions appears due to 
nonequilibrium effectsu. 

Situation is very different in the magnetic field applied 
along the layer direction, which generates the Joseph- 
son vortex lattice. Recently, dynamic properties of this 
lattice in BSCCQrhave been subject of extensive exper- 
imental researchclBBLj. At high fields (above 0.5 tesla 
for BSCGO) the Josephson vortices homogeneously fill 
all layersafl (dense Josephson vortex lattice, see Fig. [l]). 
Josephson vortex arrays in neighboring layers have strong 
inductive interactionEl, which strongly influences static 
and dynamic properties of the lattice. Due to this in- 
teraction the magnetic properties of the stack are very 
different from the magnetic properties of an isolated junc- 
tion. A transport current flowing across the layers exerts 
a Lorentz force on the Josephson vortices along the layer 
direction. An important parameter is the critical cur- 
rent above which the lattice starts to move producing 
a finite voltage. This current is determined either by 
bulk pinning or by interaction with the boundaries. In 
this paper we consider the case of a homogeneous junc- 
tions and neglect bulk pinning. The simplest and most 
known case is a single small Josephson junction with- 
out inhomogeneities, where the field dependence of the 
critical current is given by the Fraunhofer dependence, 
7 C ($) = I c0 \ sin(7r$/$o)|/(7r < i ) / < i 5 o) I with $ being the 
magnetic flux through the junction. Observation of this 
dependence has been considered as an important confir- 
mation of the dc Josephson effect. At present, the Fraun- 
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FIG. 1: Left part: Dense Josephson vortex lattice in a junc- 
tion stack. Arrows show schematically distribution of the 
Josephson and in-plane currents. Points mark centers of the 
Josephson vortices. Right part: Schematic field-size diagram 
of a Josephson junction stack. Length and field scales are set 
by the Josephson length Aj and the period of layered struc- 
ture s. This paper is focused on the long-stack/dense lattice 
regime in the upper right part of the diagram. 



hofer field dependence is routinely used as an indicator 
of the junction homogeneity. The same dependence is 
also expected for the junction stack-with the lateral size 
smaller than the Josephson lengthtil. In a single long 
junction the field dependence of the critical current has 
rather complicated-structure due to multiple coexistent 
states of the latticeO. 

In this paper we consider the critical current and lat- 
tice structure near boundaries at high magnetic fields 
in the regime of dense lattice in the stack of Josephson 
junctions. In this regime the Josephson coupling can be 
treated as a small perturbation, which allows to develop a 
full analytical description of this system. Another impor- 
tant simplification is that at high fields one can neglect 
the current self-field and neglect the difference between 
the magnetic field inside the junction and the external 
applied magnetic field. 
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We find that the boundary induces an alternating de- 
formation of the lattice. Averaging out the rapid phase 
oscillations, we derive that the lattice deformation obeys 
the Sine-Gordon equation. The deformation decays in- 
side superconductor at the typical length Is, which is 
larger than the Josephson length and increases propor- 
tional to the magnetic field, Is oc B. Structure of the 
surface deformation and the total current flowing along 
the surface is uniquely determined by the lattice position 
far away from the boundary. The total current flowing 
through the stack is given by the sum of two independent 
surface currents flowing at the edges of the sample. Due 
to the commensurability effects, the maximum current 
through the stack has oscillating field dependence, which 
resembles the Fraunhofer dependence: it has strong oscil- 
lations and overall 1/ B dependence. However, the period 
of these oscillations corresponds to adding one flux quan- 
tum per two junctions. 

We also calculate the flux-flow voltage for the Joseph- 
son vortex lattice slowly moving through a finite stack. 
Due to the interaction with the boundaries this voltage 
also has an oscillating contribution with period one flux 
quantum per two junctions. Recently, .-such oscillations 
have been observed in BSCGQ mesasLl, and also were 
seen in numerical simulations^. We find that the os- 
cillating voltage normalized to the bare flux-flow volt- 
age has universal dependence on the magnetic flux per 
one junction and the current density normalized to the 
Josephson current density. At high currents the relative 
amplitude of voltage oscillations decreases with current 
and field as 1/J 2 H 2 . We also show that in the case of 
dominating in-plane dissipation, typical for BSCCO, the 



absolute amplitude of voltage oscillations weakly depends 
on magnetic field in a wide field range, in agreement with 
experiment!! 

The new length scale Ib increasing with the magnetic 
field also sets the new field scale Bl, at which this length 
becomes of the order of the junction length L. Above 
this field the system crosses over to the regime of small 
junction and the field dependence of the critical cur- 
rent crosses over to the usual Fraunhofer dependence. 
Regimes of different behavior for a finite stack in mag- 
netic field are summarized in the phase diagram shown 
in the right part of Fig. [j]. 

The paper is organized as follows. In Section || we 
write general relations describing the phase distribution 
for layered superconductor in the high magnetic field ap- 
plied along the layers. In Section III we consider lattice 
structure near the surface, calculate the surface energy 
and surface current. In Section IV we calculate the os- 



cillating field dependence of the critical current in a fi- 
nite stack. In Section |y| we consider the field oscillations 
of the flux-flow voltage due to the interaction with the 
boundaries for slow lattice motion. 

II. GENERAL RELATIONS 



Consider a stack of the Josephson junctions with 
length L > Aj e 7s, where 7 is the anisotropy pa- 
rameter, in high magnetic field B ^> B cr = ^o/2tt^s 2 
applied along the direction of the layers (see Fig. [j]) . At 
high field the screening effects can be neglected and the 
phase distribution (p n (y) is described by the energy (per 
junction) 



dy 



dy 



2nsB 

- Ej COS [ lfn+1 — <Pn ^ V 



(1) 



where J is the in-plane phase stiffness, Ej is the Joseph- 
son energy per unit area, and w is the stack size in the 
field direction. To facilitate analysis, we introduce di- 
mensionless parameters 



= y/Aj, L = L/Xj, ej = 



£.1 



EjXjw 

The reduced energy ej is given by 
L 

1 \ - 



h = 



2n-fs 2 B 
$0 



(2) 



In a stable state the phase obeys the following equation, 
which expresses the current conservation, 



d 2 (f r , 

du 2 



+sin(<£ n+ i- ip n — hu) - sm(ip n - tp n -i — hu) = Q 



with the boundary conditions (vanishing of the in-plane 
current at the boundaries) 



dlfin 

du 



(3) 



u=0,L 



In principle, a stable state can support a finite Josephson 
current flowing through the stack 



3 = jj^jw / du sin (<p n +i - <p n ~ hu ) ■ ( 4 ) 
Jo 

In a large-size sample this current is concentrated near 
the edges. In the following section we will derive a simple 
closed expression for the edge current. 
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III. LATTICE STRUCTURE NEAR 
BOUNDARIES. SURFACE ENERGY AND 
SURFACE CURRENT. 

Far from the boundaries (exact criterion will be estab- 
lished below) the ground state configuration corresponds 
to the rhombic lattice (Fig. |l]) and the phase distribution 
can be easily calculated-using expansion with respect to 
the Josephson currentEM 



n(n + 1) 
not H 7r - 



h 2 



sin (hu — a — nn) (5) 



The bulk part of the energy is degenerate with respect 
to the phase shift a. Change of a corresponds to trans- 



lational displacement of the lattice. This degeneracy is 
eliminated by the interaction with the boundaries. To 
study this interaction, we introduce the new phase vari- 
able 4>n{ll) 



(p n = na + ir 



n(n + 1) 



and impose the condition that <p n contains only oscillat- 
ing contribution far away from the boundaries, <j) n {u) — > 
— (— l) n (2/h 2 ) sin (hu — a). This new phase <fi n (u) obeys 
the following equation 



J 



du 2 



bn+i — 4>n — hu + a + 7rn) + sin ( 

I 



hu 



(6) 



We assume that the boundary does not change the 
alternating nature of the vortex lattice and will look for 
solution of this equation in the form 

4> n {u) = (-i) n 4>{u). 

For the phase 4>(u) we obtain the following equation! 
rf 2 ( 



du 2 



2 cos 2<f> sin (— hu + a) — 



(7) 



with the boundary condition d(f>/du\ u= o t L = 0. The en- 
ergy and total Josephson current in terms of 4>(u) are 
given by 



ej(a)— / du 
Jo 



~ ( ~r~ ) ^ sin (2(f)) sin (—hu + a) 
2 \ du ' 



j(a) = —jjXj / du sin(20) cos (— hu + a) 
Jo 

. , dej 
= JJAJ^— 
oa 



, (8) 



(9) 



Note that the energy and current have the symmetry 
a — > a + ir, 4>(u) — > ~4>(u) (change of a by ir is equivalent 
to vertical displacement of the lattice by one junction). 
Therefore, the energy is 7r-periodic function of a, ej(a + 
tt) = ej(a). 

We now focus on the lattice structure near the bound- 
ary u = 0. In the limit of high field h ^S> 1 one can sep- 
arate and average out the rapidly oscillating part of the 
phase. This technique has been used in Ref. ^| to study 
melting of the Josephson vortex lattice. We split 4> into 
the smooth (v) and rapidly changing (<f>) components, 



where the smooth component describes deformation of 
the lattice induced by the boundary and satisfies condi- 
tions dv/du <C v and v — * 0, at large u. The local lattice 
compression is given by (dv/du) /h. The oscillating part 
by definition obeys the following equation 



du 2 



+ 2 cos (2v) sin (—hu + a) = 0. 



(10) 



Neglecting a weak coordinate dependence of v in compar- 
ison with the rapidly oscillating sine function, we obtain 
an approximate solution of the last equation 



— cos (2v) sin i 
h z 



-hu 



(11) 



In the limit /i> 1 the oscillating phase <j> can be treated 
as a small perturbation. Subtracting Eq. ( jToj ) from Eq. 
(0) we obtain equation for v in the first order with respect 
to 6 



d v 

——z — 40 sin (2v) sin (—hu + a) 
du 2 - 







Substituting expression for <fi from Eq. ( pi] ) and averag- 
ing out rapid oscillations, we finally obtain Sine-Gordon 
equation for the smooth lattice deformation v(u), 



h 2 d 2 v 
Tdu 2 



sin (4u) = 0. 



(12) 



Using Eq. (11), we also obtain the boundary condition 
for v(u), 



dv 
du 



du 



2 cos 2vq cos a 



(13) 



with = ^(0)- As follows from Eq. ( |l2] ) the boundary 
deformation decays at the typical length lh = h/2^/2 or, 
in the real units, 



l B = -fs 



ttjs 2 B 



(14) 



This important length scale appears as an interplay be- 
tween the shear and compression stiffnesses of the lat- 
tice. A finite system is in the large-size limit if L > 21b- 
This condition always brakes at sufficiently large field, 
B > B L , 



B L = 



L $ 
7s 2777s 2 



(15) 



Averaging out the rapid oscillations, we also derive the 
surface energy e s (a) and the surface current j s {<^) (hi 
units of jjXjw) in terms of v(u) 

1 



e s (a) ~ — sin (2v ) cos (a) 
h 



j-OO 




(dv\ 


/ du 


) 




'0 




\du ) 



1 — cos 4v 
2~h 2 



■ sin (2vq) sin (a) 



Equation ( |l2| ) has the well-known soliton solution 

tanu = tanup ex P (—2V2u/hj (16) 



Employing the boundary condition (pL3|) , we derive rela- 
tion between the boundary deformation vq and the bulk 
phase shift a 



tan(2«o) = — v2cosa 



(17) 



From this equation we can see that the maximum de- 
formation vq = =p arctan(\/2)/2 ~ =f0.478 occurs at 
a = 0,7r and the deformation vanishes at a — n/2. 
Equations (|l6| ) and (^) determine the surface deforma- 
tion v(u) for the given bulk phase shift a (i.e., lattice 
position). From the surface deformation we can restore 
all other surface properties of the lattice. Fig. || shows 
distribution of the Josephson current near the boundary 
in the two neighboring junctions for two magnetic fields, 
h = 3 and 5, and different values of the a. To illustrates 
the surface deformation, we also show the vortex cen- 
ter positions extrapolated from the bulk region assuming 
undeformed lattice. 

From the obtained solution for v(u) we derive the 
closed analytical expressions for the surface energy and 
the surface Josephson current 



e a (a) = 
js{a) = 



1 



(l - V2 + cos 2a) 



V2h 

1 sin 2a 
y/2h V2 + cos2a 



(18) 
(19) 



h=3 



h=5 




FIG. 2: Distribution of the Josephson current near the bound- 
ary it = in the two neighboring junctions for two magnetic 
fields, h — 3 and 5, and different values of the bulk phase shift 
a. Small black points mark the positions of the vortex cen- 
ters extrapolated from the bulk region assuming undeformed 
lattice. The real vortex positions deviate from these marks 
due to the surface deformation of the lattice. 



The a-dependence of the surface current is plotted in 
the left panel of Fig. ||. The maximum surface current, 
Js.max, is achieved at cos 2a = —2 + V3 and in real units 
is given by 



3s 



V3 



J J 



2-ksB' 



(20) 



Two possible directions of this current correspond to the 
Lorentz force on Josephson vortices directed from the 
boundary and towards the boundary, promoting either 
entrance or exit of vortices. In contrast to the regime of 
dilute lattice at B < B cr , in which the maximum values 
of the entrance and exit current are very different (see, 
e.g., Ref. |l4|), in the dense lattice they are approximately 
the same. It is interesting to note that, up to numeri- 
cal factor, the expression for the entrance current in the 
dilute lattice regime coincides with Eq. (p0|). 



IV. CRITICAL CURRENT OF A FINITE STACK 

Consider now the total current flowing through a finite 
stack. For the given phase shift a (lattice position) the 
total current, Jl(«) = Jl(o, hL, h), is the sum of the two 
surface contributions coming from the boundaries u = 
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FIG. 3: Left Panel: Plot of the function hj s = -L "j" 2 " , 

J J V2 v2+cos 2a ' 

which gives dependence of the surface current on the phase 
shift a (see Eq. (|l9|)). Right Panel: Dependence of function 
/ = hj c , which determines the shape of the critical current 
oscillations in stack of long junctions, on the reduced flux 
through one junction \ = 27r<J?/<J>o in Eq. (0). 
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FIG. 4: Thick solid line shows dependence of the critical cur- 
rent for the stack of long Josephson junctions on the magnetic 
flux per junction. For comparison we also show by dashed line 
the Fraunhofer dependence |sin(7r<J> / 3>o) | /(7r$/$o) which de- 
scribes the field dependence of the critical current in a single 
small Josephson junctions. 



Jl{q)- 



sin 2a 



sa\2{hL — a) 



y/2h W2+cos2a y/2+cos2(hL-a) 



(21) 



Introducing notation \ = hL = 27r$/$o, where 4> = BsL 
is the magnetic flux through one junction, we obtain that 
the maximum current «/ c (Xj h) = fix)/h is given by 



Jc{x, h) = max [J L (a, x, h)} 



(22) 



Numerically obtained dependence of f — hJ c vs x is 
plotted in the right panel of Fig. [| In the real units, the 
field dependence of the maximum Josephson current in 
the field range $o/2ti"7s 2 < B < L$ /7T7 2 s 3 is given by 



2-ksLB 



f 



2-ksLB 



(23) 



where Ij = jjLw is the maximum Josephson current 
through the stack at zero field, and oscillating func- 
tion f(x) is plotted in the right panel of Fig. To 
facilitate comparison with experiment we also obtain 
an interpolation formula for fix) f° r < \ < tt/2, 
f(x) ~ 0.128 + 0.888 cos(x) + 0.021 cos(3x)- The de- 
rived field dependence of the critical current somewhat 
resembles the well-known Fraunhofer dependence in a 
single small-size Josephson junction: it has overall 1/B 
dependence with large oscillations. However, it has a 
very different physical origin and also has several quali- 
tative differences. The most important difference is that 
the period of oscillations corresponds to adding one flux 
quantum per two junctions, i.e., it is two times smaller 
than for the Fraunhofer oscillations. This kind of oscil- 
lations have been recently observed in the flux-flow volt- 
age of BSCCO mesasQ. Also, contrary to the Fraunhofer 
dependence, the points where the magnetic flux inside 
the junction equals to integer flux quanta, correspond to 



the local maxima of the critical current. Dependence of 
the critical current on the magnetic flux per junction is 
plotted in Fig. [|. For comparison, the Fraunhofer depen- 
dence is also shown. The dependence (|2^) holds until two 
edges give independent contributions to the total current. 
This condition breaks when the magnetic field exceeds 
B L = L$ /(27T7 2 s 3 ). At higher fields the field depen- 
dence crosses over to the usual Fraunhofer dependence. 
Recently, this crossover has been studied by numerical 
simulationsEj. 



V. OSCILLATING FLUX-FLOW VOLTAGE 

When the external current exceeds the critical current, 
the lattice starts to move. At slow motion the surface de- 
formation has time to adjust to the current lattice posi- 
tion. In this case the surface energy produces a periodic 
potential for the moving lattice and one can use static 
results to predict the I-V dependencies. Time variation 
of the lattice phase shift obeys equation 



da . 
l/ff ~dt + ^' = ' 



(24) 



where the total surface current Jl(cv) = Jl(q;, hL, h) is 
given by Eq. ( |2l| ) (for brevity we again skip in equations 
its dependence on the magnetic field and size) and the 
viscosity coefficient, Vff, is related to the flux- flow resis- 
tance of the stack, Rff, 



v ff 



N<I>o 
2ncRff 



where N is the number of junctions in the stack. Volt- 
age drop per one junction U is related to da/dt by the 
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FIG. 5: Oscillating dependencies of the reduced voltage on the 
magnetic flux per one junction for different current densities. 
The plot represents a universal behavior for junction stacks 
with different sizes. 



Josephson relation 

_ $ da 

2?7C dt 

Solution of Eq. ( p4[ ) is given by the implicit relation 

f u ff da ' =f 
Jo J-J L (a>) ' 

from which we obtain the average phase change rate 

v Sf da 



da 
~dt 



1 



Jo J — Jl{o) 
and the flux-flow voltage 

da 



U 
Wf 



J 



7T J J - J L {a) 



i -i 



(25) 



with Uff = RffJ being the bare flux- flow voltage with- 
out periodic potential. Because the surface current 
Jl{o) = Jh{oi,hL,K) oscillates with the magnetic field, 
this flux-flow voltage will also experience similar field 
oscillations. Such oscillations have been recently ob- 
served by Ooi et. alu Using scaling property of the cur- 
rent Ji(a, hL, h) = LF(a, hL) we can see from Eq. ( p5| ) 
that the reduced flux- flow voltage U/Uff depends only 
on two parameters: the current density normalized to 
the Josephson current density, i = J/L = j/jj, and the 
magnetic flux through one junction, $ = BsL. This al- 
lows for a unified description of the voltage oscillations 
in junction stacks with different sizes. Oscillating depen- 
dencies of U/Uff on the magnetic flux $ for different 
current densities are shown in Fig. ||. 

At large currents one can expand this result with re- 
spect to surface current and obtain a small correction to 



the flux-flow voltage 
U 

Uff 



1 - 



g(Q) + g(x) 

TTj 2 h 2 



(26) 



where the oscillating part is described by dimensionless 
function 



g(x) 



sin 2a sin 2 (a — x) 



da 



o V2 + cos 2a^2 + cos 2(x - a) tt 
0.263 cos 2x + 0.0046 cos4x- 



Therefore, the absolute amplitude of voltage oscillations 
is given by 



U.A,) 2 (-. 



(27) 



Independently from the dissipation mechanism, the rel- 
ative amplitude of oscillations 5U/RffJ decreases with 
field as 1/B 2 . In the regime of dominating in-plane dis- 
sipation^ypical for BSCCO, the flux-flow resistance is 
given byt£l 



Rff = R c 



B 2 



B 2 + B 2 



B n 



$ 



\[2-k^ 2 s 2 



where a a b and a c are the components of quasiparticle 
conductivity. Using this relation one can rewrite Eq. j2^ ) 
as 



SU = 



RabJ] 



1 



2 J 1 + B 2 /B 2 



(28) 



with J M = c^ /8n 2 X 2 b , R ab = Ns/Sa a b, and S = wL 
is the junction area. For BSCCO, typically, B a ~ 10T. 
In the region B < B a the oscillation amplitude weakly 
depends on the magnetic field, in agreement with exper- 
iment (see Fig. 2 in Ref. It interesting to note that 
in this regime it is mainly determined by the in-plane 
parameters of superconductor. 



VI. SUMMARY AND ACKNOWLEDGEMENTS 

In conclusions, we found that the edge current of the 
dense Josephson lattice is uniquely determined by the 
magnetic field and position of the lattice in the bulk. 
Near the surface the lattice has alternating deformation, 
which decays inside superconductor at the typical length 
which is proportional to the magnetic field. Due to the 
rhombic lattice structure, both the critical current and 
the flux-flow voltage at small velocities have oscillating 
field dependencies with the period of one flux quantum 
per two junctions up to the size-dependent magnetic field. 
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